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Study guide 
Introduction 


1 Differential equations and direct integration 
1.1 Differential equations and solutions 
1.2 Direct integration 


2 Solution by separation of variables 
2.1 Further use for the Chain Rule 
2.2 Separation of variables 


3 Modelling growth and decay 
3.1 Radioactive decay 
3.2 Continuous model for population growth 


4 Differential equations with the computer 
5 A look ahead 


Summary of Chapter C3 
Learning outcomes 


Summary of Block C 
Solutions to Activities 
Solutions to Exercises 


Index 


There are five sections in this chapter. They are 
intended to be studied consecutively in four study 
sessions. Section 4 requires the use of the 
computer and Computer Book C. 


The pattern of study for each session might be as 
follows. 


Study session 1: Section 1. 
Study session 2: Section 2. 
Study session 3: Section 3. 
Study session 4: Section 4 and Section 5. 


Fach session requires two to three hours. 

Sections 2 and 3 may take longer to study than the 
other sections. Section 5 is short and will not be 
assessed. 
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Introduction — 


In many modelling situations, the ultimate mathematical aim is to obtain 
an explicit relationship between two variables. This involves expressing the 
dependent variable, which we denote here by the symbol y, as a function F 
of an independent variable, denoted by x. However, it may not be possible 
to write down an equation of the form y = F(z) straight away. This desired 
outcome may be achievable only by solving some other type of equation. 


Often the equation that we have to solve involves not only the variables x 
and y themselves but also derivatives of y with respect to x. You saw 
examples of such situations in Chapter C2, Section 3, in the context of 
modelling motion. Equations that include derivatives, such as 

2 
ay =. 62 ¥ + 9x*y = sin(2z), 
are called differential equations. Any function y = F(x) which satisfies 
such an equation is called a solution of the differential equation. For 
example, the function y = e*” is a solution of the first differential equation 
displayed above. 


The order of a differential equation is the order of the highest derivative 
that appears in it. So the first equation displayed above is an example of a 
first-order differential equation, since it contains only the first derivative of 
the dependent variable with respect to the independent variable. A 
second-order differential equation, such as the second equation displayed 
above, contains a second derivative, and maybe a first derivative as well, 
but no higher derivatives. 


This chapter concentrates exclusively on first-order differential equations, 
and only on particular types of these. Section 1 deals with differential 
equations that can be written in the form 


where f is a given function. Such equations can be solved by direct 
integration of the function f. Section 2 shows how differential equations 
that can be written in the form 


= = f(x)g(y) 


may be solved, using the method of separation of variables. This method 
again depends on being able to perform integrations. 


Section 3 looks more closely at a special type of differential equation 
soluble by separation of variables, namely dy/dxz = Ky, where K is a 
constant. This type of equation can be used to model decay and growth 
processes, such as radioactive decay or the growth of a population. 


Section 4 shows how the computer can be used to plot direction fields, 
which give a visual impression of the overall behaviour of solutions of a 
differential equation. This is followed by a demonstration of how solutions 
of differential equations may be obtained both numerically and graphically, 
using EKuler’s method. These approaches apply more widely than those 
considered earlier in the chapter. 


Finally, Section 5 indicates briefly some directions in which the theory of : 
differential equations can be further extended. 


In these examples the 
independent variable was t, 
denoting time, while the 
dependent variable was 

v (velocity) or s (position). 


In this chapter we normally 
use Leibniz notation for 
derivatives. This is common 
in the context of differential 
equations. 


Higher-order differential 
equations are defined 
similarly. 


1 Differential equations and direct integration - 


Solutions expressed in the 
form of ‘a more general 
equation relating x and y’ are 
discussed in Section 2. 


See Chapter C2, Section 3. 


These rules were introduced 
in Chapter Cl, Section 4. 


In Subsection 1.1 we look at first-order differential equations and their 
solutions. In Subsection 1.2 the method of direct integration is used to 
solve a simple type of differential equation. 


1.1 Differential equations and solutions 


We start by drawing together some of the terminology referred to already 
in the Introduction. 


Definitions 


© A differential equation is an equation that relates an 
independent variable, x say, a dependent variable, y say, and one 
or more derivatives of y with respect to x. 


The order of a differential equation is the order of the highest 
derivative that appears in the equation. A first-order 
differential equation involves the first derivative, dy/dx, and no 
higher derivatives. 


A solution of a differential equation is a function y = F(z) (ora 
more general equation relating x and y) for which the differential 
equation is satisfied. 


As you saw in Chapter C2, a differential equation may feature symbols for 
the independent and dependent variables other than x and y, and you will 
see further examples of this later in the current chapter, when a modelling 
context makes it appropriate. While considering the general theory of 
differential equations, however, we shall stick to the use of x and y. 


A first-order differential equation features the first derivative, dy/dx. The 
differential equation may also involve separate occurrences of x or y, or 
both. Thus each of the following is a first-order differential equation: 
ies a St ee = ep 
dx dx dx 3y?+1 
You will shortly see methods which lead to the solutions of all these 
differential equations. First, however, you are asked to verify in a number 
of cases that a function which is claimed to be the solution of a differential 
equation does in fact satisfy that equation. This type of exercise is useful 
in checking that no error has been made in the working leading to a 
solution. It is similar to the practice of checking whether a calculated 
indefinite integral is correct by differentiating the answer obtained. 


In addition, such checking will provide you with some further practice in 
differentiation, using Leibniz notation, including application of the 
Product, Quotient and Composite Rules. The following example indicates 
what is required. 
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Example 1.1 Verifying that a given function isa solution | 


In each case, show that the given function satisfies the given differential 
equation. 


d 
fet ee se + 22: sige ee, 
dx 
x dy 1 — x’ 
b = ~——: SS = 
(b) y 1+ 2?’ dx (1+27) 
d 3+ 2x 
(c) y= ae"; x = — (x > 0) 
d 
(d) y=In(1 +2”); 2s = 2re7¥ 
Solution 


In each case, we differentiate the given function to find an expression 
equivalent to the left-hand side of the differential equation, and then show 
that the same expression arises from the right-hand side of this equation. 
This verifies that the function is a solution of the differential equation. 


(a) The differential equation is dy/dx = 3x* + 2. The derivative of the 
given function, y = x° + 22, is 

ds 

— 22) = 327 +2 

q, (t+ 2a) = 3x" + 2, 

which is the same expression as the right-hand side of the differential 

equation, as required. 


(b) The differential equation is dy/dx = (1 — x*)/(1+ 27)’. In order to 
differentiate the given function, y = x/(1+ 27), we apply the Quotient 
Rule. The function can be written as y = u/v, where 


u=x and v=142’, 


for which 
du dv 
—— eS 
is 1 — - x 
Hence we have 


dy 1 ( du =) 
aso" % de dx 
(1+ x*)(1) — (w)(22) 
(1+ 27)? 

LG Papa? so 2pF85 1 mg" 

Se tt)” 
which is the same expression as the right-hand side of the differential 
equation, as required. 


(c) The differential equation is dy/dx = (3+ 2x)y/zx. In order to 
differentiate the given function, y = x°e*”, we apply the Product Rule. 
The function can be written as y = uv, where 


u=a2> and v=e”*, 

for which 
du dv 
—=32? and — =2e”. 
dx ae = 


We describe each function F' 
by writing y = F(z). 


See Chapter Cl, — 
Subsection 4.2. 


See Chapter Cl, 
Subsection 4.1. 


The restriction x > 0 in the 
statement of part (c) ensures 
that division by z is valid 
here. Another suitable 
restriction is x < 0. 


See Chapter Cl, 
Subsection 4.3. 


Recall that expa is an 
alternative (and in this case 
convenient) way of writing e°, 
and that exp is the inverse 
function of In. 
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Hence we have 
dy du dv 
de dx | dx 
= (3a°)(e*”) + (x*)(2e*) = (3 + 22)z*e?*. 
Now we substitute y = r%e?” into the right-hand side of the differential 
equation, to obtain 


3 22 
(3 + 2x)y = (3 + 2r)x°e = (34 Dene. 
x z 
This is the same expression as was obtained above for the derivative, 
as required. 
(d) The differential equation is dy/dx = 2xe~¥. In order to differentiate 
the given function, y = In(1 + 27), we apply the Composite (Chain) 
Rule. The function can be written as 


y=lInu, where u=1+ 2’, 
for which 

ay. Tt du 

a an ee 

a dx 


Hence we have 


dy _ dy du 
dx dudz 
1 22 
aS ) 1+ 2? 


Now we substitute y = In(1 + 2?) into the right-hand side of the — 
differential equation, to obtain 
2re-¥ = 2x exp(—1n(1+ 2”)) 
= 2x exp(In((1 + x*)~*)) 


PA 
= IJ¢(1 +27) ' = 
( ) oa 
This is the same expression as was obtained above for the derivative, 
as required. 


Here are some similar cases for you to try. 
Activity 1.1 Verifying that a given function is a solution = 


In each case, show that the given function satisfies the given differential 
equation. 


dy i) 

(a) y = 5/2; as oe (x > 0) 
3 dy 3 3 
(b) y = sin(e**); =. 3e°” cos(e**) 
£ 
d 

(c) y =e cost Ze): = =y(1—2tan(22)) (—f7 < x < $n) 

l-¢ dy 
d ram : ee eee ‘ 2 
Cees Se 


Solutions are given on page 43. 
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In Example 1.1(a), we showed that the function y = x° + 2z is a solution 
of the differential equation dy/dx = 3x” + 2 because 


aC 4 Qz) = 32° + 2. 
x 


However, this is not the only solution of the differential equation. Since 


=(2° + 2e-4+--c) = 327 + 2, 

where c is any constant, any function of the form y = x° + 2x + c is also a 
solution of the differential equation dy/dx = 3x7 + 2. We say that 

y = x? + 2x +c, where c is an arbitrary constant, is the general solution of 
this differential equation, since it describes the whole infinite family of 
possible solutions. By contrast, the function obtained when a particular 
value is chosen for c, as with the original function y = x° + 2x (where 

c = 0), is called a particular solution of the differential equation. 


You may recall that an arbitrary constant also arises when finding the 


indefinite integral of a function and, as you will see shortly (or may For a discussion of integrals 
already have noticed), there is a very close link between integration and and arbitrary constants, see 
the process of solving differential equations. Chapter C2, Section 1. 


The general solution of a first-order differential equation will usually 
feature an arbitrary constant. This does not always occur in the form of a 
‘te’ term tagged on at the end of the solution formula, but that will be 
the pattern for the differential equations solved in this section. 


General and particular solutions 


© The general solution of a differential equation is the set of all 
possible solutions of the equation. It usually involves one or more 
arbitrary constants. 


A particular solution of a differential equation is a single 
solution of the equation, which consists of a relationship between 
the dependent and independent variables that contains no 
arbitrary constant. 


In many cases, where a differential equation arises in a mathematical 
model, the value of the dependent variable is known for one value of the 
independent variable. For example, in Chapter C2 the motion of a See Chapter C2, Example 3.1. 
parachutist was modelled by a differential equation for dv/dt, where v is 
velocity and ¢ is time. In addition, a value of v was specified at time t = 0. 
This permitted an appropriate value to be chosen for the arbitrary 
constant that arose from integration, and hence led to a single function (a 
particular solution) which described the velocity of the parachutist. 


Similar considerations apply more generally. A first-order differential 
equation has a general solution and, in order for a particular solution to be 
obtained from this, a further condition on the solution is needed. This 
often takes the form of an initial condition, in which a value of the 
dependent variable y is specified at a given value of x. When taken 
together, the differential equation and initial condition are called an 
initial-value problem. 
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The second notation stems 
from the common practice of 
writing 


y = y(z), 


in which the symbol y does 
‘double duty’, as a variable on 
the left-hand side and as a 
function on the right. This is 
convenient because it reduces 
the number of separate 
symbols that might otherwise 
be required. Nor does it cause 
any ambiguity, since the 
context determines which of 
the two uses is intended. 


This is revision of what you 
saw in Chapter C2, 

Example 3.1, Activity 3.1 and 
Subsection 3.2. 


This can be seen by adapting 
the solution to 
Activity 1.1(a). 
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Initial-value problem 


© An initial condition associated with a first-order differential 
equation requires that the dependent variable y takes a specified 
value, 6 say, when the independent variable x has a given value, a 
say. This is often written as 


y=bwhenz=a, oras y(a)=b. 


The numbers a and 0 are called initial values for x and y, 
respectively. 


The combination of a first-order differential equation and an 
initial condition is called an initial-value problem. The 
solution of an initial-value problem is a particular solution of the 
differential equation which also satisfies the initial condition. 


The following example and activity concern finding a particular solution 
from the general solution of a differential equation, using an initial 
condition. 


Example 1.2 Applying an initial condition 


As you have seen, the general solution of the differential equation 
dy/dx = 327 +2 is y= 2° +2x+¢, where c is an arbitrary constant. Find 
the particular solution that satisfies the initial condition 


Y= 5 whee v= 2. Ses. 9) = 9. 
Solution 


We substitute the values x = 2 and y = 9 (simultaneously) into the 
expression y = 2° + 2xr+c. This gives 


Ga 42x otc Dts 
so c= —3. The required particular solution is therefore 


y He ae Ss 


Here is a similar problem for you to try. 


Activity 1.2 Applying an initial condition 


The general solution of the differential equation 

dy 5 

Se 1 

tbe fae ( ) 
is y = 5,/x +c, where c is an arbitrary constant. Find the particular 
solution that satisfies the initial condition 

y =ob whee bj thet: is; Ol) = 


A solution is given on page 43. 


1.2 Direct integration 


In the remainder of this section we concentrate on differential equations 
which can be written in the form 


—= f(z). (1.1) 


As you saw in Chapter C2, differentiation is ‘undone’ by integration, and 
so the general solution of the differential equation (1.1) is given by the 
indefinite integral of the function f(x); that is, 


y= [#@ dx. 


Example 1.3 Finding the general and a particular solution 
(a) Find the general solution of the differential equation 
dy = 
oe =—wW. 
(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 2 when x = 1. 


Solution 


(a) The general solution of the differential equation is given by 


y= f ode. 


The general solution is therefore 
y=3r +e, 
where c is an arbitrary constant. 


(It is always possible, and advisable, to check a general solution by 
differentiating and substituting back into the differential equation. 
Here y = 32° + c, so dy/dx = 2”, as required.) 


(b) Putting x = 1 and y = 2 into the general solution, we obtain 
2= 2(1)’ +e. 
Hence c = e, and the required particular solution is 


SS SS 
y= {fF =: 


The method used in Example 1.3 is called the direct integration method 
for solving differential equations. It can be applied to any differential 
equation, such as dy/dx = x’, for which the derivative of the dependent 
variable is equal to a known function of the independent variable. 


See Chapter C2, Section 1. 
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Direct integration 


© The general solution of the differential equation 


dy 
—— f(z) 
is the indefinite integral 
Recall that F(x) is an = [#@ dx = F(x) +¢, 
integral of f(z) if 
F where F(x) is any integral of f(z) and c is an arbitrary constant. 
aa (x)) = f(a). Any initial condition 
Hence y = F(z) isa y =:.b. when. = a,.. thatis,.. #(0)= 6, 


particular solution of the 


differential equation enables a value for the arbitrary constant c to be found. The 


corresponding particular solution satisfies both the differential 
equation and the initial condition. 


The method can also be applied, of course, if symbols other than x and y 
are used for the independent and dependent variables. Any differential 
equation of the form 


d(dependent variable) 3 
ditndeeden’ Gara dependent bl 
d(independent variable) f (independent variable) 


We say ‘in principle’ because can in principle be solved using direct integration. For example, the 


in practice it may not be differential equation 
possible to find explicitly an d 

integral of f. However, you -— = Af +3 

will not meet such cases in dt 

this text. has the general solution 


q= | (4t+3)dt= 2 + 3t+e 


where c is an arbitrary constant. 


Activity 1.3 Can direct integration be applied? 


You are not asked to solve To which of the following differential equations can the direct integration 
these differential equations. method be applied? 
d dx dx 
(a) = = EC08 2 (b) a = ROS (c) FH = boost 
dy dx dx 
d) —2=ycos e) — = f) —=zcost 
(d) 7 = ycosy (2) fae ycosy — (f) = = xcos 


Solutions are given on page 43. 
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Activity 1.4 Applying direct integration | 


(a) Find the general solution of the differential equation 
dy 22 
Pe =—€E . 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 2 when x = 0. 


Solutions are given on page 44. 


In Activity 1.4 you saw that the general solution of the differential equation 
dyjds =e is y = se" +c, where ¢ is an arbitrary constant. Some graphs 
corresponding to this family of solutions are shown in Figure 1.1. Each of 
these solution curves has the property that, at any point (x,y) on it, the 
gradient of the curve is dy/dx = e**. All the solution curves in Figure 1.1 
are parallel, in the sense that the vertical separation of any two curves is a 
constant, independent of x. This is a property of all differential equations 
of the type dy/dx = f(x). However, as you will see in the rest of the 


chapter, it is not a general property of first-order differential equations. 


Activity 1.5 Applying direct integration again 


Solve the initial-value problem 


=sint, x2£=1whent=0. 


dt 
A solution is given on page 44. 


Activity 1.6 Mass of a burning rocket 


A toy rocket consists of a casing and fuel. Initially the total mass of the 
rocket is 200 grams. The fuel is ignited at time t = 0 and then burns in 
such a way that the rocket’s total mass m grams at time t seconds satisfies 
the differential equation 


(t > 0). 
When the fuel is exhausted, the residual mass of the rocket casing is 
100 grams. 


(a) Find the general solution of the given differential equation. 


(b) Find the particular solution of the differential equation that satisfies 
the initial condition m = 200 when t = 0. 


(c) Hence find the time at which the fuel is exhausted; that is, find the 
time t when m = 100, for the particular solution obtained in part (b). 


Solutions are given on page 44. 


You can check the general 
solution by differentiating it 
and then substituting back 
into the differential equation. 


Figure 1.1 Graphs of 
y = 5e7" +c, forc = —3, 0 


3 
and 5 


In this activity the 
independent variable is t and 
the dependent variable is z. 


The derivative dm/dt 
represents the rate of change 
of mass. This quantity is 
negative because the mass 
decreases as the fuel is burnt. 
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Summary of Section 1 


This section has introduced: 


© first-order differential equations, and the solutions of such equations 
(general and particular); 


© the process of obtaining a particular solution from a general solution, 
using an initial condition; 


© the method of direct integration, to solve differential equations of the 


form dy/dx = f(z). 


Exercises for Section 1 


Exercise 1.1 
Show that the function y = —1/z (x > 0) is a solution of the differential 
equation dy/dx = y’. 
Exercise 1.2 
(a) Find the general solution of the differential equation 
dy 
—_— = > 0). 
(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 5 when z = 4. 


(c) Hence find the value of y when z = 25. 


Exercise 1.3 


Solve each of the following initial-value problems. 


d 
(a) = = cos(2r), wu =—2 when x = in 
os 1 


Exercise 1.4 


A large spherical snowball melts in such a way that, at any instant, the 
rate of decrease of its volume V is proportional to its surface area A, with 
constant of proportionality k > 0. Initially the snowball has a radius of 
100cm. After 2 days the snowball has lost half its initial volume. 


Using the formulas V = ¢ar° and A = 4rr? for the volume and surface 


area of a sphere, in terms of its radius r, it can be shown that the change 
in radius r (in cm) of the snowball with time t (in days since melting 
started) is modelled by the differential equation 


dr 
You are asked to establish = —k (r>0). 
this differential equation in t 
Exercise 2.4. (a) Find the general solution of this differential equation. 


(b) Find the particular solution of the differential equation that satisfies 
the given initial condition. 


(c) Use the condition that the snowball has lost half its initial volume 
after 2 days to find the value of the constant of proportionality k. 


(d) How long will it take for the snowball to disappear completely? 
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ylutio separation of variables 


In Subsection 1.2 you saw that the general solution of the differential 
equation 


dy 
Sees f(z) 
is given by the indefinite integral 


y= [#@ dx. 


In this section we develop a method for finding solutions of a further type 
of differential equation, namely, those of the form 


— = f(x)g(y), 


where f and g are known functions. This method, called separation of 
variables, is described in Subsection 2.2. As a preliminary, Subsection 2.1 
shows how the Composite (Chain) Rule for differentiation may usefully be 
employed even when one of the component functions in the composite is 
not known explicitly. 


2.1 Further use for the Chain Rule 


In Section 1 you saw the solutions of several differential equations. In each 
case, the solution was specified in the form y = F(x), for some function F’. 
When solving a differential equation, it is sometimes not possible to reach 
such an explicit form for the solution function. For example, suppose that 
it has been deduced, from a differential equation for dy/dz, that the 
variables x and y are related by the equation 


y ty=2+9. (24) 
It is not straightforward to solve this equation for y and hence to reach the 
form y = F(x) for the solution. However, equation (2.1) does describe a The possibility of a solution 
solution, even if only implicitly. For example, a corresponding solution being ‘a more general 


curve can be plotted, and the value of dy/dzx at a point (z,y) on this curve equation relating x and y’ was 


is the slope of the curve at that point. allowed for in the definition of 
a solution of a differential 


More generally, suppose that the solution of a differential equation is equation in the box on page 6. 
obtained in the form 


H(y) = F(2), (2.2) 


where H and F are known functions, and H(y) # y. (Equation (2.1) is of 
this form.) This is called an implicit solution of the differential equation. 
By contrast, a solution written in the form y = F(x), where F is a known 
function, is called the explicit solution of the differential equation. In 
seeking to solve a differential equation, it is the explicit solution which we 
seek. However, if this is impossible to find, then we must settle for an 
implicit solution. 


In Subsection 1.1 you saw how to verify that a given function is an explicit See Example 1.1 and 
solution of a specified differential equation. In order to check that a Activity 1.1. 
calculated implicit solution satisfies the differential equation from which it 

was obtained, the Chain Rule is employed. This is possible because, in 

equation (2.2), the y that appears on the left-hand side is (implicitly) a 


1 


function of xz. (All the implicit solutions we deal with have the form of 
equation (2.2) and are such that y is a function of x.) If we write 


Here z is an additional se iiiy) = Fi). 
variable, introduced only to 


iastik the cieminedt at this then two expressions can be deduced for dz/dz as follows. Since z = F(z), 


point. we have 
dz 
— = F(z). 
= (x) 
Since also z = H(y), where y is a function of x, we have, by the Chain 
Rule, 
dz - “aera dy 
— = —— =A’ (y)—. 
dx dydzx ) dx 
(This step is called implicit differentiation of H(y) with respect to z.) 
Equating these two expressions for dz/dz, we see that 
d 
if H(y) = F(z), then H'(y)— = FU (2.3) 
* 
The next example indicates how this result can be put to use. 
Example 2.1 © Verifying a solution given implicitly | 
(a) Show that equation (2.1), 
y ty=2°+9, 
is an implicit solution of the differential equation 
dy 2x 
dx 3y?+1 
Note that (1,2) does lie on (b) What is the slope of the corresponding solution curve at the 
the given curve, since point (1,2)? 
2°+2=10=17+49. Sesletiean 


(a) We have H(y) = F(x), where H(y) = y?+y and F(x) = x? + 9. It 
follows from result (2.3) that 


d d 
= se —(s" ¥ oF; 


d 3 


dy 
which leads to 
dy 
= 
Dividing through by 3y? + 1 gives 
dy 22 
dx 3y? +17 
as required. 
(b) At (1,2), the slope of the curve is 
dy a 2 
ae Foes 


(3y* + 1) Qn. 


Here are similar problems for you to try. 


| | 16 


SECTION 2. SOLUTION BY SEPARATION OF VARIABLES 


Activity 2.1 Verifying a solution given implicitly 


(a) Show that the equation 
y+siny=2+e"—1 
is an implicit solution of the differential equation 
d 1 +e” 
eer rege a 
(b) What is the slope of the corresponding solution curve at the 
point (0,0)? 


Solutions are given on page 44. 


Activity 2.2. Another verification 


In Chapter C1 you saw that 


—(tan 2) = sec’ 2. 
7, tan) 
Use this result, and the trigonometric identity sec? 6 = 1 + tan’ 0, to show 
that the equation 
tany=2z 
is an implicit solution of the differential equation 
ee 
dx 1+2? 
A solution is given on page 44. 


1 1 
(5% a E < aad. 


Comment 


In fact, it is possible in this case to write the solution explicitly. For if 
tany = x and —i0 << ¥< ST, then y = arctan. Hence the result of this 
activity shows that | 


“ (arctan 2) 

—(arctan x) = ————. 

dx 1+ 2? 

It can be shown in a similar way that 
1 


VJ/1 — x? 


Gy (aresin £) = (-l<2< 1). 


The previous example and activities indicate the usefulness of implicit 
differentiation using the Chain Rule, as described by result (2.3). In the 
next subsection, you will see that result (2.3) enables us to explain a 
solution process for a wider class of differential equations than those 
considered in Subsection 1.2. 


See Chapter Cl, 

Activity 4.4(a). The same 
differentiation appears in 
Chapter C2, Activity 2.4(a). 


This trigonometric identity 
was introduced in 
Chapter C2, Subsection 2.1. 


This derivative formula is 
valid for all values of zx. 
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The phrase ‘in implicit form’ 
is used to describe a general 
or particular solution that is 
an implicit solution. The 
phrase ‘in explicit form’ is 
applied similarly. 


This demonstrates that, when 
both sides of an equation are 
integrated, we need introduce 
only one constant of 
integration. We shall follow 
this practice from here on. 


2.2 Separation of variables 


We now derive a method for solving differential equations of the form 


where f and g are known functions. The starting point for explaining this 
method is the result, obtained in Subsection 2.1, that 
dy 
then A'(y)— = F(z). 
(y) 5 (x) 
Integration is the reverse process to differentiation. Hence, faced with a 
differential equation of the form 


if Aly) = fe), (2.3) 


where f and A are known functions, a solution will be given by the 
equation H(y) = F(a), provided that 


H"(y)=h(y) and F(x) = f(a). 
It follows that 
if n(y) => f(z); “thet [rw) iy = [f@ dx. (2.4) 


The latter equation gives the general solution of the differential equation, 
in implicit form. 


The next example shows how this result may be applied. 


Example 2.2 Finding a general solution in implicit form 


(a) Find, in implicit form, the general solution of the differential equation 
dy 
a 
7 dz 
(b) Find the corresponding particular solution that satisfies the initial 
condition y = 3 when zx = 2. 


327. 


Solution 


(a) Comparing the differential equation with result (2.4), we have 
h(y) = 2y and f(x) = 3x”. Hence the general solution of the 
differential equation is given by 


| 2ydy = [30° ax; that is, y?+a=2° +5, 


where a and 6 are arbitrary constants. These two constants can be 
combined to form the single arbitrary constant c= b— a. Thus the 
general solution, in implicit form, is 

y == +¢, 
where c is an arbitrary constant. 


(b) To apply the initial condition, put x = 2 and y = 3 into the general 
solution. ‘his gives 


7 =2 +< 
so c= 1. Hence the required particular solution is 


y= a +1. 


Here is a similar problem for you to try. 


Activity 2.3 Finding a general solution — 


(a) Find, in implicit form, the general solution of the differential equation 
<~—=-1 (y>0). 


(b) Find the corresponding explicit form of this general solution. 


Solutions are given on page 45. 


The approach demonstrated in Example 2.2 and Activity 2.3 can be 
extended to solve any differential equation of the form 


a = S (99). v=) 


A differential equation of this type is said to be separable, and the 
corresponding method of solution is called separation of variables. ‘To 
solve equation (2.5), we first rewrite it so that there is a function of x alone 
on the right-hand side. This is done by dividing both sides by g(y), to give 


= f(a). 


g(y) dx 


Now we can apply result (2.4), with h(y) = 1/g(y). This shows that the 
general solution of equation (2.5), in implicit form, is 


/ ay t= / Son ay: (2.6) 


To summarise, the method of separation of variables is as follows. 


Separation of variables 
The method applies to differential equations of the form 


= F(z)a(y). 


1. Divide both sides by g(y), to obtain 
L dy _ 


g(y) dx = f(x). 


2. Integrate both sides with respect to x. From result (2.4), the 


outcome is 


[a5u- [ #@ dx. (2.6) 


Carry out the two integrations, introducing one arbitrary 
constant, to obtain the general solution in implicit form. If 
possible, manipulate the resulting equation to make y the subject, 
thus expressing the general solution in explicit form. 


(A general solution, in implicit or explicit form, can be checked by 
substituting back into the differential equation. ) 


ARATION OF VARIABLES 


We assume here that any 
value y for which g(y) = 0 is 
excluded from consideration. 


At this stage the variables are 
separated, with a function of 
y alone on the left-hand side 
and a function of x alone on 
the right-hand side. 


19° 


The combination of a 
differential equation and an 
initial condition is called an 
wnitral-value problem — see 
page 10. 


Some differential equations of 
the type given by 

equation (2.7) will be 
considered in Section 3. 
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Any initial condition provides a value for the arbitrary constant in step 3. 
The corresponding particular solution satisfies both the differential 
equation and the initial condition. 


You might like to note that the equation in step 2 can be written down 
quickly from the differential equation obtained in step 1 by regarding 
dy/dx temporarily as an ordinary quotient and ‘multiplying through by dz’ 
before adding the two integral signs. This shortcut is a convenient aid in 
recalling this step of the method. However, remember that the derivative 
dy/dx is not dy divided by dx. We have not assigned a meaning to either 
dx or dy on their own; they make sense only when combined in a derivative 
such as dy/dx, or when preceded by an integral sign as in [... dz. 


This method can be applied also when independent and dependent 
variables other than x and y are used, provided that the derivative is equal 
to a product of a function of the independent variable and a function of 
the dependent variable. 


Activity 2.4 Can separation of variables be applied? 


To which of the following differential equations can the separation of 
variables method be applied; that is, which of them are of the form 


ee = f(independent variable)g(dependent variable) 
d(independent variable) 


for some choice of the functions f and g? For each of the equations below 
that is of this form, identify the functions f and g. 
dy dx dp _ 
(a) qe 7 ECOSY (b) 7H = teost (c) 7 PcosP 
d d 
@) Faery (©) Ga=vtey  (f) F =coslay) 


Solutions are given on page 45. 


Parts (b) and (c) of Activity 2.4 illustrate that both the equations 


oY _ §(2) (1.1) 


= S600) al 


are special cases of the equation dy/dx = f(x)g(y), with g(y) = 1 in the 
first case and f(x) = 1 in the second. As you saw in Subsection 1.2, the 
general solution of equation (1.1) can be found by direct integration, to 
which separation of variables reduces when g(y) = 1. 


Example 2.3 Applying separation of variables 


(a) Find the general solution of the differential equation 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 1 when z = 0. 


Solution 


(a) 


The differential equation is of the form of equation (2.5) with 


f(x) =cosz and g(y) =e ¥. Dividing both sides of the differential Other choices for f(x) and 
equation by e ¥ gives g(y) are possible, for example, 
Fiz} = 2c and 
dy 
e’— = cose. g(y) = se ¥, but all such 
choices lead to the same 
From this we obtain, using result (2.4), general solution. 
Since 1/e-¥ = eY, dividing by 
vas BRE. 
J dy = [ cosa da. e 4 is the same as 


we . 
On performing the integrations and including a single arbitrary — 
constant c, this becomes 

e’ =sinzx+c, 
which is an implicit form of the general solution. 


To make y the subject of this equation, recall that In(e”) = y. Hence, 
on taking the natural logarithm of each side of the equation, we obtain 


y = In(sinz +c), 


where c is an arbitrary constant. This is the explicit form of the 
general solution. 


Putting xz = 0 and y = 1 into the explicit form of the general solution, 
we have 


1 = In(sin0 +c) = Ine. 
Hence c = e’ =e, and the required particular solution is Recall that In(e) = 1. 


y = In(e+sinz). 


Comment 


‘; 


The solution could now be checked, though the details are not given 
here. | 


In part (b) the explicit form of the general solution was used in finding 
the required particular solution. This particular solution can also be 
obtained by using the implicit general solution. 


Putting x = 0 and y = 1 into 
e’ = sinzx+c, 
we obtain 
e=sin0+c, 
so c=e. Hence the particular solution in implicit form is 
e’ =e+sing, 
and taking the natural logarithm of each side gives 
y = In(e+sinz), 
as before. 


Unless you are directed otherwise, you should find particular solutions 
from the explicit form of the general solution. 


Recall that sec x = 1/(cosz). 
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Here are some similar problems for you to try. 


Activity 2.5 _ Applying separation of variables 


(a) Find the general solution of the differential equation 


dx y 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 1 when x = 0. 


Solutions are given on page 45. 


Activity 2.6 Applying separation of variables again 


(a) Find the general solution of the differential equation 


d 
o = tsecx (—sn < x < $n). 


dt 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition xz = 0 when t = 0. 


Solutions are given on page 46. 


A possible snag in applying the separation of variables method was 
mentioned in a margin note on page 19. We can divide both sides of the 
equation 

d 

i = f(x)gy) 
by g(y) only if g(y) 4 0, since division by zero is an undefined operation. 
It is therefore necessary to exclude from consideration any value of y, the 
dependent variable, for which g(y) = 0. This may be done by imposing a 
suitable condition on y from the outset, as we did in Activity 2.5 by 
assuming that y > 0. 


As well as having to impose restrictions on the dependent variable in order 
to apply the separation of variables method, we must be aware of other 
pitfalls that can arise. The solution obtained may be valid only for a 
certain range of values of the independent variable. For example, in 
Activity 2.5 you found that the required particular solution of the given 
differential equation dy/dx = —x/y, assuming that y > 0, was 


y= V1—2?. 


The expression \/1 — x? is a real number only for values of x such that 
—1<a< 1. Since x = 1 and x = —1 both give y = 0, which is excluded by 
the restriction y > 0 in the question, the domain of this particular solution 
is the interval (—1,1). Similar considerations apply to the corresponding 
general solution in Activity 2.5, which was 


y = Vb— x?. 


This formula provides positive values for y only when the arbitrary 
constant b is positive, and then only when —b!/? < x < b/?. So the 
interval (—b'/?, b!/2), where b > 0, is the domain of this solution. 


SECTION 2. SOLUTION BY SEPARATION OF VARIABLES” 


These examples show that consideration of the appropriate domain for a 
solution can become a bit complicated. This should not be a cause for 
your concern, as you will not be expected to derive the appropriate domain 
for a solution in the manner indicated above. However, it is as well to be 
aware that in practice the choice of domain may require some thought. 
Deriving a formula y = F(x) (or H(y) = F(«)) for a solution does not 
guarantee that the formula obtained makes sense for all real numbers 2, 
nor that the function F actually satisfies the differential equation without 
some restriction being imposed on its domain. Also, as has been pointed 
out several times, it is wise to check a proposed solution by substituting it 
into the differential equation. 


The final activity of this section investigates a model for the motion of a 
space probe launched from Earth. The modelling theme is continued in the 
next section, where each model features a differential equation that can be 
solved by separation of variables. 


Activity 2.7 Motion of a space probe 


A motorless space probe is launched vertically upwards from the Earth’s 
surface. When the probe is at a distance r from the centre of the Earth, its 
outward velocity vu satisfies the differential equation 


dv — — (eS #, The restriction r > FR arises 
dr oF because the differential 
where R is the radius of the Earth and k is a positive constant. All the equation is a suitable model 
variables and constants are measured in the appropriate SI units. only when the space probe is 

above the surface of the 


(a) Motion of the probe outwards from the Earth corresponds to v > 0. 


Find the general solution of the differential equation in this case. 


(b) The initial condition corresponding to a launch velocity of vp is 
v(R) = vp; thatis, v=vr whenr = R. 


Show that the corresponding particular solution for the motion 
outwards from Earth is | 


1 1 


(c) Show that at very large distances from the Earth, the velocity of the Do not spend too long on 


space probe tends towards the limiting value part (c) if you find it difficult. 
2k 
2 
Up —- —. 
Pave 


Hence, by considering the cases vg > \/2k/R and vp < V2k/R, 


explain the significance of the particular launch velocity 


Solutions are given on page 46. 
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Summary of Section 2 


This section has introduced: 


© the contrast between the explicit solution of a differential equation, in 
the form y = F(x), and an implicit solution, in the form H(y) = F(z); 

© implicit differentiation of a function H(y) with respect to x, where y is 
a function of x whose rule is not known explicitly; 


© the method of separation of variables, to solve differential equations of 
the form dy/dx = f(x)g(y). 


Exercises for Section 2 


Exercise 2.1 
(a) Show that the equation 


y+Iny = (1—2z)e** 
is an implicit solution of the differential equation 


dy y(1—6z)e** 
de pakke tT 


(b) What is the slope of the corresponding solution curve at the 
point (0,1)? 


Exercise 2.2 


(a) Find the general solution of the differential equation 


du 
2 4 
u—=Z2. 
dx 
(b) Find the particular solution of this differential equation that satisfies 
the initial condition wu = 1 when zx = 1. 


Exercise 2.3 


Solve in explicit form each of the following initial-value problems. 
d 
(a) = 9° iy > F,. y = > when «=D 
es 


d 
(b) — = te (230), = 2 wer tab 


Exercise 2.4 


In this exercise you are asked to establish the differential equation 


dr 
—-=-—k > 0 
bo ee 


which was solved in Exercise 1.4. The equation models the change in 
radius r (in cm) of a snowball with time t (in days since melting started). 


The snowball melts in such a way that, at any instant, the rate of decrease 
of its volume V = ¢rr® is proportional to its surface area A = 4rr?, with 
constant of proportionality k > 0. Use these formulas for V and A, and the 
Chain Rule, to establish the differential equation. 


SECTION 2. SOLUTION BY SEPARATION OF VARIABLES 


Exercise 2.5 


This exercise concerns a mathematical model for the flow of water from a — One application of this model 
tank, of constant horizontal cross-section, through a small hole near its is to the water emptying from 
base. According to the model, the height h of the water surface above the _ vertically-sided bath via the 
hole is related to the time t since the hole was uncovered. This relationship plug-hole. 

is described by the differential equation 


dh 
dt — =f (h = 0), 


where k is a positive constant. (The value of k depends on both the area of 
the hole and the horizontal cross-sectional area of the water column in the 
tank.) 


(a) Find the general solution of this differential equation. 


(b) Find the particular solution for which h = ho at t = 0. (Here ho is the 
height of the water column above the hole at the moment when the 
hole is uncovered.) 


(c) Show that, according to this model, the tank empties to the level of 
the hole after a time 2hq/7/k. 
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3 Moding pom oet'anay 


In this section we investigate mathematical models that are based on 
first-order differential equations of the form 


se Ee 
oF Y; 


where K is a constant. This type of equation is soluble by separation of 
variables. The resulting behaviour is that of exponential growth or decay. 


Subsection 3.1 concerns radioactive decay, while Subsection 3.2 returns to 
the topic of population change. 


3.1 Radioactive decay 


As you may know, the nuclei of radioactive substances disintegrate 
spontaneously, at a rate which varies according to the substance concerned. 
For example, a gram of pure uranium is composed of an enormous number 
(about 2.5 x 10°") of uranium atoms. Each of these atoms consists of a 
central nucleus surrounded by a cloud of 92 electrons. There are three 
different kinds of naturally-occurring uranium nuclei, called uranium—234, 
uranium—235 and uranium—238, but uranium as mined is more than 99% 
uranium—238. (The numbers 234, 235 and 238 give the atomic mass of 
these three kinds of uranium.) From time to time one of the uranium—238 


In Chapter B3, nuclei disintegrates, emitting an electrically charged particle, called an 
Subsection 4.1, a particle was alpha particle, which can be detected using a Geiger counter. As a result 
defined as a material object of these disintegrations, the uranium slowly decays into another substance, 
whose size and internal namely, thorium—234. The amount of uranium contained in a lump of 
structure may be neglected. matter therefore decreases steadily with time. 

The usage here is taken from 

physics, where ‘particle’ Each nucleus of uranium—238 has a known probability of decay in a given 
means an object which is time interval. However, because of the enormous number of nuclei present 
finite but very small. in even one gram of uranium—238, it is sensible to use a continuous model 


to describe this situation, measuring the amount of uranium in grams, 
rather than using a discrete model, where we would measure the number of 
uranium nuclei. The same applies to the decay of other radioactive 
substances. 


As each nucleus has a constant probability of decay in a given time 
interval, we assume for the continuous model that the rate of decay of a 
radioactive substance is proportional to the amount of that substance 
which is present. Consequently, the change in mass m of the radioactive 
substance that is present at time t can be modelled by the first-order 
differential equation 


dm 

—=-km (m>0O), Sl 

= (m > 0) (3.1) 
If time ¢ is measured in days, where k is a positive constant, called the decay constant. The negative 
then the unit of k is day~?. sign in this differential equation reflects that we have a decaying process, 


where the mass m is a decreasing function of time, so that dm/dt is 
negative. ‘The condition m > 0 corresponds to the fact that decay can 
occur only when some of the radioactive substance is present. Physically, 
we would expect equation (3.1) to be a good model for the situation, 
except where there are only a small number of radioactive nuclei, for which 
a discrete model would be more appropriate. 
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To find a particular solution of equation (3.1), we require an initial 
condition. If there is an amount mo of the radioactive substance present at 
time t = 0, then the initial condition is 


m=m whent=0; that is, m(0) = mo. (32) 


The next activity asks you to solve the initial-value problem given by 
equations (3.1) and (3.2). 


Activity Cpe Solving the radioactive decay problem 


(a) Use separation of variables to find the general solution of 
equation (3.1), giving your answer in explicit form. 


(b) Find the corresponding particular solution which satisfies the initial 
condition (3.2). 


Solutions are given on page 46. 


In the above activity you showed that the initial-value problem that 
models radioactive decay, 


d | 
—- = km (m > 0), Mm = ™Mp when t = 0, 
has solution 
m =e. (3.3) 
a Figure 3.1 Graph of 
The graph of this solution, shown in Figure 3.1, illustrates that the meme k > 0 


solution has the characteristics we would expect of radioactive decay. The 
mass of the substance decreases, initially quite quickly and then more 
slowly, until eventually there is practically none of the substance left. 


Different radioactive substances have different values for the decay 
constant k. The larger the value of k, the faster the substance decays. 
This is illustrated in Figure 3.2. 


The exponential decrease of a radioactive substance has an interesting 
consequence in terms of the factor by which the amount of the substance 
declines over any fixed time interval. If we denote the mass of the 
radioactive substance at time t by m(t) then, according to equation (3.3), 
we have 


= —kt 
ees Figure 3.2 Graph of 
After a further time interval of duration T, the mass of the radioactive m = moe“ for several 
substance has become values of k 
mi P= Winer d= Mee ee = MES. Here we use the rule 


It follows that et TY — ered, 


m(t + T) =F 
a (3.4) 


This equation says that, at the end of a time interval of duration T’, the 
proportion of the original amount of a radioactive substance remaining is 
e-*? regardless of the time t at the start of the interval. 


According to equation (3.4), the time which it takes for the mass of 
radioactive substance to diminish to half its original amount is given by 


ee m(t + T) xg *F. 


: m/(t) 


Recall that 
In 5 = In(2~*) = —In2. 


You do not need to be able to 
follow the detail of the 
following outline description 
in order to attempt 

Activity 3.3. 


In terms of numbers of atoms, 
this proportion is 1 to 10!?. 


This discovery was made in 
the 1950s by the American 
chemist Willard F. Libby. He 
was awarded the Nobel Prize 
for Chemistry in 1960 for his 
pioneering work in 
radiocarbon dating. 
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Taking the logarithm of each side, we obtain 


In 5 Ri; that is. Ki = in 72. 


Hence the mass of radioactive substance decays to half of its original 
amount in the time 

pm? 

k; 

where k is the decay constant. This time T is called the half-life of the 
substance. Whereas the speed of decay of a radioactive substance could be 
described in terms of the decay constant k, it is more usual to use the 
half-life T. Radioactive substances have a very wide range of half-lives, 
from millions of years down to tiny fractions of a second. 


(3.5) 


Activity 3.2. Using the half-life 


Silicon—31 has a half-life of 2.62 hours. Find the value of the decay 
constant k for silicon—31. Hence calculate the proportion of a sample of 
silicon—31 which will still be present after six hours. Give each value to 
three significant figures. 


A solution is given on page 47. 


Carbon dating 


One useful application of radioactivity is to the dating of archaeological 
finds of biological origin. Carbon occurs naturally in three forms: 
carbon—-12, carbon—13 and carbon—14. Chemically they are 
indistinguishable, but carbon—14 is radioactive, with a half-life of 
approximately 5570 years, and decays continually (into nitrogen). This 
decay is balanced by the production of carbon—14 in the atmosphere, and 
it may be assumed that the proportion of carbon—14 in naturally-occurring 
carbon is constant over time. This same known proportion exists in all 
living tissue. 


When an organism dies and is buried, the carbon—14 is no longer 


replenished and the proportion of carbon—14 in the carbon decreases with 


time. The quantity of carbon—14 in a sample of dead tissue can be 
estimated from a measurement of the radioactivity, and the total quantity 
of carbon in the sample can be estimated from chemical analysis. These 
values, the known half-life of carbon—14 and the known proportion of 
carbon—14 in naturally-occurring carbon in living tissue form the basis for 
a method of estimating when the organism died. The method is known as 
carbon dating. In particular, an estimate of the proportion of the 
original amount of carbon—14 still present in the sample can be found. 


Activity 3.3 Using carbon dating 


The bones of an animal are found in an archaeological dig. Analysis 
produces the estimate that 85% of the original amount of carbon-—14 is still 
present in the bones. Taking the half-life of carbon—14 to be 5570 years, 
find the approximate age of these bones to three significant figures. 


A solution is given on page 47. 
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Finding the decay constant from data 


You have seen that the amount m of a radioactive substance present in a 
lump of matter at time t is modelled by an initial-value problem whose 
solution is 


Mm = Moe,". (3.3) 


The decay constant k cannot be measured directly, but it can be deduced 
by measuring experimentally how the amount of the radioactive substance 
changes with time. The data in the table below were obtained by Ernest 
Rutherford early in the twentieth century, when he was studying 
radioactive decay. He observed a sample of a radioactive compound of 
thorium over a period of several months, taking measurements from time 
to time, and deducing the proportions m/mg of the original mass of the 
sample which remained. 


O41 hie Ui G46 8.33 “0.22 ©.13 0:06 


According to the model, the amount of radioactive substance present, as a 
proportion of its initial mass, is given by m/mp = e~™'. Taking the 
logarithm of each side of this equation, we obtain 


= (= | scree 
Mo 


Hence a plot of In(m/mpo) against t, for the data above, should result in a 
straight line through the origin with slope —k. After adding to the table 
above a row for the values of In(m/m,) to two decimal places, we obtain 
the table below. 


t (days) 10 15 20 26 33 48 66 89 
mmo 671i <—-858- 051 —0.43.-500993 =: OO = O85 0.08 
iia i ae 0 67 084-111 =15] —204 281 


In Figure 3.3 a graph of In(m/mo) against t has been plotted. It shows 
both the data points and a straight line through the origin which fits them 
well. ‘The line is drawn through the origin because, when t = 0, we have 
In(m/mo) = In(mo/mo) = Inl1 = 0. 


In(m/mo) 20 AO 60 80 100 + 


Figure 3.3 Graph of In(m/mo) against t for the data in the table above 


mg is the original amount of 
the substance. 


Ernest Rutherford 
(1871-1937) was a New 
Zealand-born physicist of 
great distinction. He was the 
first to show that 
radioactivity is caused by the 
spontaneous decay of heavy 
atoms into slightly lighter 
ones, for which he was 
awarded the Nobel Prize for 
Chemistry in 1908. He also 
propounded the nuclear 
model for the atom, and 
became the first person to 
split the atom. 


A straight line through the 
origin in the (xz, y)-plane, with 
slope —k, has the equation 

y = —ka. Here we have t in 
place of x and In(m/mo) in 
place of y. 


In Block D you will see an 
objective method for 
calculating the line which 
best fits a set of data points. 
However, it is often sufficient 
to use one’s eye. This is a 
rather subjective approach, 
but a good rule of thumb is to 
draw the line so as to 
minimise the distances of the 
data points from the line, 
while ensuring that there are 
about the same number of 
data points on either side of 
the line. This rule of thumb is 
employed in Figure 3.3. 
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From the graph you could 
read the second coordinate as 
slightly different from —3.2. 
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Figure 3.3 shows that these data points are reasonably well-fitted by a 
straight line through the origin with negative slope. This validates use of 
the exponential model to describe radioactive decay. The straight line 
passes through the origin and the point (100, —3.2), so its slope is 
approximately —3.2/100 = —0.032. Hence the decay constant for 
radioactive decay of the thorium compound is estimated to be 


k= 0082 day. 
From equation (3.5), the corresponding half-life is 


T In2 ry In 2 
k 0.032 


A number of remarks about the above estimation, the decay constant and 


the half-life are called for. 


1. It is assumed in the above method for estimating the decay constant 
that hand plotting of the data and the graph are used. 


= 22, Gays {to 2 s.I.). 


2. The values of In(m/m,) were given to two decimal places because hand 
plotting on a normal-sized piece of graph paper can represent this level 
of accuracy at most. | 

3. The first coordinate of the point chosen in order to calculate the slope 
of the line drawn by eye should be near the high end of the t data 
range, and ideally be arithmetically convenient. The first coordinate of 
the point (100, —3.2) satisfies both requirements. 

4. The value of JT’ was given to two significant figures, the same accuracy 
as for the value of k, estimated from the graph. Slightly different 
estimates for k, obtained from other lines ‘fitted by eye’, can affect the 
value of T. For example, the point (100, —3.25) leads to k = 0.0325 
and T' = 21.3 to three significant figures. 


You should bear the above points in mind when attempting the following 
activity, and in similar situations. 


Activity 3.4 Finding the half-life of uranium-—239 


The table below shows data for the radioactive decay of a sample of 
uranium—239. 


t (minutes) 


0.73 0.56 0.42 0.29 0.24 


(a) Calculate to two decimal places the entries for the third row of this 
table. | 

(b) Plot,.on graph paper, a graph of In(m/mo) against t for these data. 
Draw ‘by eye’ a line through the origin which you think best fits the 
data points. 


(c) Using your straight line from part (b), estimate values for the decay 
constant k and the half-life T for uranium—239. 


Solutions are given on page 47. 


SECTION 3 MODELLING GROWTH AND DECAY — 


Comment 


As was said before, fitting a straight line to data points ‘by eye’ is a 
subjective process. This means that your values of the decay constant and 
half-life may not be exactly the same as those given in the solution. 
However, they should not be very different! 


3.2 Continuous model for population growth 


Earlier in the course you saw mathematical models, based on recurrence Discrete models for 
relations, to describe the way in which population size alters with time. population change were 
These models all had the property that they predicted population size only developed in Chapter B1, 
at discrete intervals of time, typically at annual intervals. The aim here is Sections 2~4. 

to develop a model which describes a population at any instant of time, on 

a continuously changing basis. In practice, of course, populations are made 

up of individuals, and so the size of a population can take only integer 

values and does not change continuously. However, if the population size is 

large, then its behaviour may be approximated satisfactorily by a 

continuous function, and continuous models for populations are therefore 

of value. 


If P represents the size of a population at time ¢, then the instantaneous 
rate of change of the population size is given by the derivative dP/dt. A 
simple modelling assumption is that the rate of change of the population 
size is proportional to the current population size. We then have the 
differential equation 
dP 
ie: ARF {eo 
where K is a constant, called the-proportionate growth rate. The 
condition P > 0 corresponds to the fact that population change can occur 
only for a positive population. 


The value of K may be positive (birth rate exceeds death rate), negative 
(birth rate less than death rate) or zero (birth and death rates equal). If 
the initial population size (at time t = 0) is Po, then we have the 
initial-value problem 


dP 7 
This problem has the same mathematical form as that given by The only difference is that —k 
equations (3.1) and (3.2) for the radioactive decay model. Hence the had to be negative (k > 0), 


solution to problem (3.6) can be written down immediately, by comparison whereas K (in place of —k) 
with the previous case whose solution is given by equation (3.3). The can take any real value. 
solution to problem (3.6) is therefore 


- 
re Pye**. (3.7) 


Figure 3.4 shows some typical graphs of this solution. The case K > 0 
corresponds to population growth, K = 0 gives a static population size and 
K <0 corresponds to population decline. The behaviour for K < 0 is 
mathematically identical to that examined earlier for radioactive decay 
(where K = —k), so we shall now look at examples for which kK > 0. The 
next activity shows how equation (3.7) can be put to use. 


Figure 3.4 Graph of 
P = Poe** for several 
values of Kk 


si 


See Chapter B1, 
Subsection 2.2. 


See Activity 2.2 of 
Chapter B1. See also 
Figure 2.6 of that chapter. 
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Activity 3.5 Exponential growth of a human population — 


The population size of a country is currently 5 million, and its 
proportionate growth rate is 3% per year, so that K = 0.03 year~!. 


(a) Use equation (3.7) to predict what the population size of the country 
will be after ten years and after a hundred years. 


(b) After how long will the population size reach 50 million? 


In each part, give your answers to three significant figures. 


Solutions are given on page 48. 


For obvious reasons, this model for population change is called the 
exponential model. You may recall that in Chapter B1 this name was 
given to the first discrete model for population change. Which of the two 
models is being used should be clear from the context, but if there is any 
possibility of confusion, then we distinguish between them by referring to 
either the discrete exponential model or the continuous exponential model. 


As already noted, the continuous exponential model for populations has 
the same mathematical form as that of the radioactive decay model. It 
follows that, in place of equation (3.4), we now have 


PPD orang 
abies sak (3.8) 


This equation says that, in a time interval of duration T, the population 
size changes by a factor e*7, regardless of the time t at the start of the 
interval. If A < 0, then the population is in decline and has a half-life T, 
given by T = —(In2)/K, just as for radioactive decay. If K > 0, then the 
population size increases. Here we can define the doubling time T for the 
population, which is the time that it takes for the population to double in 
size from any starting value. 


Activity 3.6 Doubling time — 


(a) Suppose that the growth of a population is described by a (continuous) 
exponential model, with proportionate growth rate K > 0. Determine, 
in terms of AK, the doubling time for the population. 


(b) Find the doubling time, to three significant figures, for the population 
described in Activity 3.5. 


Solutions are given on page 48. 


Finding the parameters of the exponential model 


As with the decay constant for radioactive decay, the proportionate growth 
rate for a population may be deduced from data on population sizes. In 
the case of a human population, this will be census data. 


For example, the table below shows census data for the USA for the years 
1790-1890. We analysed these data to an extent in Chapter B1, using the 
discrete exponential model. There we estimated the (discrete) 
proportionate growth rate by using two data points, namely, the 
population sizes in 1790 and 1890. 
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This has the advantage of relative simplicity, but a more accurate answer 
can be expected if all the available data are taken into account. We shall 
now do so in the context of the continuous exponential model, to estimate 
a value for the proportionate growth rate K. 


Year 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890 The population is given in 
Pome a 5 Oe 2S EE 8S 314 398-052 62.9 Millions, to the nearest 
100 000. 
We start from the equation 
Putte, (3.7) 
and take the logarithm of each side. This gives 
ln P = In(Poe*"*) . Recall that In(e”) = 2, and 
= In Py + In(e**) = In Py + Kt. that In(zy) = Inez + Iny. 


Now compare this equation with the general form of the linear function 

y = mz-+c, as follows: Here m means the constant 
slope of a linear function, and 
not the mass of a sample of 

t : t radioactive substance. 


ln PoE iat, 


| eee 2 eee eae © 


This comparison shows that, provided the exponential model is Similar comparisons were 
appropriate, we have a linear relationship between the variables In P and t. carried out for the radioactive 
Hence a plot of In P against t should result in the data points lying on or — decay model, and log-linear 
close to a straight line which crosses the (In P)-axis at In Pp and has plots were drawn from data. 


slope K. Such a plot of In P against t is called a log-linear plot. In that case, c = 0. 


Example 3.1 US population (1790-1890) 


(a) By using a log-linear plot for the data in the table above, show that 
the exponential model is appropriate for the US population for the 
years 1790-1890. Measure time ¢ in years, from t = 0 at the 1790 
census date. 


(b) Use your plot to estimate values for the proportionate growth rate Kk 
and the initial population size Pj in the exponential model. 


(c) Use this exponential model to predict the population of the USA 
in 1950. 


Solution 


(a) The table below gives the values of the time t, the population size P 
and In P, to three significant figures. 


t 0 10 20 30 AO 50 60 70 80 90 100 


r 39: 53 T2 96 129 171-232 34 2 oe 6 
in?’ * 136 1.67 “1.07 226. 256 24 314 345-306 597° 4.14 


The corresponding log-linear plot is shown in Figure 3.5. 
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The straight line was drawn 
‘by eye’. 


~The value of the parameter Po 
for the model differs from the 
actual value of the population 
at time t = 0 because the 
model value is chosen by a 
process of line-fitting that 
involves all the given data for 
the population. 


In P Ire 


: 20 40 60 80 100 # 
Figure 3.5 Plot of InP against t for the data in the table 


The data points are fitted reasonably well by the straight line drawn in 
the figure. This confirms that the exponential model is appropriate to 
describe the US population for the years 1790-1890. 


The intercept of the line on the (In P)-axis is 1.40; this is the value of 
In Pp for the model. Since In Pp = 1.40, we have 

Py =e =400 00 
As well as the point (0, 1.40), another point on the straight line is 
(100, 4.23). The slope of the straight line is therefore 

4.23 — 1.40 

100 — 0 

So our estimate for the proportionate growth rate K is 0.0283 year~', 
and for the initial population size Pp the estimate is 4.06 (million). 


= 0.0283. 


Using the estimates obtained for K and P,, the exponential model for 
the US population (1790-1890) is 


PP ae 4.06e°:0283¢ | 


The year 1950 corresponds to t = 160, so our estimate for the 
population in that year is 


P(160) = 4.06¢e°-°783*16 — 376 (to 3 s.f.). 
So our estimate for the population of the USA in 1950 is 376 million. 


Comment 


The actual US population in 1950 was only 150.7 million! This shows that 
the exponential model used here cannot be expected to give good results 
much beyond the period 1790-1890 for which the data were given. 


In the following activity, you can carry out a similar exercise for another 
population. 


The table below shows data for the population of India (in millions, to the 
nearest million) for the years 1950-2000. 


1950 1960 1970 1980 1990 2000 
398 442 555 689 851 1014 


Year 


Population 


(a) By using a log-linear plot for the data in this table, show that the 
exponential model is appropriate for the Indian population for the 
years 1950-2000. Measure time t in years, from t = 0 at the 1950 
census date. 


Use your plot to estimate values for the proportionate growth rate K 
and the initial population size Pp in the exponential model. 


(b) 


(c) Use this exponential model to predict the population of India in 2020. 


Solutions are given on page 48. 


Mathematical note 


For each of the radioactive decay and population models, we restricted the 
dependent variable (m and P, respectively) to be positive. In terms of the 
independent variable x and dependent variable y, we have derived the 
result that dy/dx = Ky (y > 0) has general solution y = Ae**, where A is 
a positive but otherwise arbitrary constant. The restriction on y arose 
from modelling considerations, though it was also relied on in the method 
of solution. 


If the restriction on y is removed, then the expression above for the general 
solution remains correct, but with no restriction on A; that is, 
y= Ae“*, 


has the general solution (3.9) 


where A is an arbitrary constant (which may be positive, negative or zero). 
This solution can be checked, as usual, by verifying that the differential 
equation is valid when y = Ae*® is substituted into it. 


Comparison with the discrete exponential model 


This subsection concludes by considering how the continuous exponential 
model for population change compares with the corresponding discrete 
model developed in Chapter B1. The discrete model was given by the 
recurrence relation 


Foy = (5+ 


where P,, is the population size at n years after a specified initial time and 
r is the (discrete) proportionate growth rate (giving the factor by which 
the population increases year on year). This recurrence relation has the 
closed-form solution 


P, =(1+1r)"Po 


(n = 0,1,2,...), 


(n= 0,1, 2, ..-). 
From either the recurrence relation or its solution, we have 


Past 
Py, 


3 47. (3.10) 


[ODELLING GROWTH AND DECAY 


The method used the fact 
that 


1 
[i dy=tny +e 
] 


for y > 0. An alternative 
approach is required to deal 
with the cases y < 0. 


The remainder of this 
subsection will not be 
assessed. 


See Chapter B1, 
Subsection 2.2. 


For example, if K = 0.03 then 
r= 0.0305, and if K = 0.2 
then r ~ 0.22. 


The discrete model does, to 
the extent that it focuses on a 
‘snapshot’ at annual intervals. 


For example, the continuous 
logistic model, which is the 
analogue of the discrete 
model in Chapter B1, 
Sections 3-4, is given by the 
differential equation 


dP FP 
—=KP| t= — 
dt ( a) 


where K and EF are constants. 
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This can be compared with the corresponding formula for the continuous 
model, which is 


PUt+T) _ xr 
= ee: (3.8) 


Now P,, and P,,,, are population sizes taken one year apart, at t = n and 
t=n-+1. Hence we can match equations (3.8) and (3.10), for integer 
values of t, by putting t = n and T = 1 into equation (3.8). We see then 
that the two models will give the same values for P(n) = P,,, provided that 


e*=1[+4+r. 


This equation relates the discrete proportionate growth rate r (measured 
over one year), with the continuous proportionate growth rate K 
(measured instantaneously). For a given population, these two 
proportionate growth rates are not quite the same. However, if the 
population size varies slowly, then the two values will be close to each 
other. 


As was pointed out, a continuous model would be inappropriate for a 
population whose size is small. Continuous models have the apparent 
advantage over discrete models that they predict the population at all 
times, and not just annually. However, for many species the birth rate is 
highest in the spring whereas the death rate is highest in the winter. The 
continuous exponential model does not take account of this seasonal 
variation. 


Furthermore, as discussed in Chapter B1, for very large populations the 
environment becomes less hospitable. Birth rates then tend to decline and 
mortality rates to rise, so in practice the proportionate growth rate 
decreases with the size of the population. There are continuous models, 
based on differential equations, which take account both of seasonal 
variations and of the decrease in proportionate growth rate with population 
size. However, there is not the space to study them in this course. 


Summary of Section 3 


This section has introduced: 


© the (continuous) exponential models for radioactive decay and for 
population change; 

© amethod for finding parameter values in these models from data, 
using a log-linear plot; 


© the form of differential equation dy/dx = Ky (occurring in each of the 
above models), with general solution y = Ae*?”. 
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Exercises for Section 3 


Exercise 3.1 


The interest on bank deposit accounts is usually compounded annually, Compounding interest 

six-monthly or quarterly. Suppose that, in order to attract customers, a n times per year, at a rate 

certain bank introduces a type of account on which interest is compounded 7” per year, means that an 

continuously at the rate of 5% per year. = credit balance cp rises 
O 


(a) Use the fact that the proportionate rate of increase of any credit a 
balance is a = 0.05 to write down a differential equation for the Co € = 4 
credit balance c (in £) of such an account at time t (in years), during 


° . . ° i - ft : 
any time interval in which no deposits or withdrawals are made. ee 
(b) Write down the general solution of this differential equation. _ 


(c) A customer places £100 in such an account and makes no further 
deposits or withdrawals in the subsequent year. Find the 
corresponding particular solution of the differential equation, and 
hence calculate the credit balance of the account after one year. 


(d) Compare the interest received by this customer after one year with 
that which would be obtained at a rate of 5% per year compounded 


(i) annually; (ii) six-monthly; (iii) quarterly. 
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4 Differential equations with the computer 


Es 


For each point (x,y) in the 
(two-dimensional) domain 
of f, a unique real value 


f(x,y) is defined. 


To study this section you will need access to your computer, together with 
Computer Book C. 


In Section 2, you saw how to solve equations of the form dy/dz = f(x)g(y) 

using the method of separation of variables. There are many first-order 

differential equations which cannot be expressed in this form. For example, 
dy 


a, 4.1 
dx Ss? (4.1) 


is one such equation. More generally, we can consider differential equations 
of the form 


dy 


rs = f(x,y). (4.2) 


Here f is a known function of two variables, x and y. In some cases such 
differential equations can be solved to obtain a formula (either implicit or 
explicit) for the solution. Such a solution is said be analytical. On the 
other hand, it is sometimes not possible to find a formula which describes 
a solution. In such cases, the behaviour of functions that satisfy a 
differential equation can be deduced instead by numerical means, using a 
computer. You will see how this is done here. 


We start by considering the graphical information which may be extracted 
from a first-order differential equation. A solution of equation (4.2) will 
have a graph which is a curve in the (z, y)-plane. For each point (z, y) that 
the curve passes through, the gradient of the curve has the corresponding 


value f(z, y). 


However, the differential equation specifies such gradient values 
throughout the domain S, in the (2, y)-plane, of the function f, and not 
just on any particular curve that passes through S. Hence the behaviour of 
the differential equation (as opposed to any individual solution of it) can 
be represented graphically if we have a means of visualising gradients at all 
points of the (2, y)-plane. 


The association of a gradient with each point in the domain S is called a 
direction field for the differential equation. The gradient at each point 
can be represented satisfactorily by a short line segment that has the 
appropriate gradient and is centred on the point concerned. It is not a 
practical proposition to draw such line segments at all points of S, but we 
can do so for a representative grid of points. If the spacing between grid 
points is not too large, then this will provide a good picture of how the 
differential equation behaves. 


The direction field for the differential equation (4.1) is represented in 
Figure 4.1. It can be seen that the gradients shown do correspond to the 
numerical values of f(x,y) = x+y. For example, the gradient is zero at 
the origin, —1 at (—2,1) and 2 at (1,1). 
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Figure 4.1 Direction field for dy/dx =x +y 


The graph of any solution to the differential equation must have a gradient 
at each point which matches that provided by the direction field. You can 
probably see that such graphs could readily be sketched on top of the 
direction field diagram, ‘following the flow’, to give a rough idea of how 
solutions behave. There are also more accurate ways of constructing such 
graphs, as you will see shortly. 


Refer to Computer Book C for the work in this section. 


You have seen how solutions to a differential equation of the type (4.2) can 
be obtained numerically. In fact, each such solution corresponds to a 
particular initial-value problem, 


Y f(x,y), yle0) = wo 
This can be solved using Euler’s method, which is described by the pair The Swiss mathematician 
of recurrence relations Leonhard Euler (1707-83) 
was hugely prolific. He 
Ent1 = In +h, Yat = Un +hf(tn,yn) (n=0,1,2,.. .). devised the computational 


method that now bears his 
name in order to calculate the 
orbit of the Moon. It was 
used by the German 
cartographer and astronomer 


In these equations, h is the step size between the successive values of x at 
which solution estimates are calculated. Each calculated value y, is an 
estimate of the corresponding ‘true solution’ y at x = 2, that is, y, 1s an 
estimate of y(x,,). The sequence of estimates depends on the choice of 


both the step size h and the overall number of steps N. Decreasing h, (Johann) Tobias Mayer 

while increasing NV to cover the same range of z-values, leads to (1723-1762) to construct 
progressively improved estimates for the solution values, and with a small —_jynar tables. In 1765, Mayer’s 
enough step size, any desired level of accuracy can be achieved. widow received £3000 from 


the British Board of 
Longitude for the 


Summa ry of Section 4 contribution which his tables 


made to the problem of 
This section showed how first-order differential equations can be visualised determining longitude, while 
graphically, using direction fields, and how solutions may be obtained Euler received £300 for his 


either numerically or graphically, using Euler’s method. a contribution to the 
| work. 
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5 A look ahead 


This section will not be 
assessed. 


You can check this solution 
by substitution. 


This is called the simple 
harmonic motion equation, 
and provides the simplest 
model of an oscillating 
system. Again, you can check 
the solution by substitution. 


For example, z might be the 
temperature at the point 
(x,y) of a flat metal plate. 


For example, the variation in 
temperature O(z,t) at 
distance xz along a thin rod at 
time ¢ is modelled by the heat 
equation, 

070 | ee 

da? At? 


where k is a constant. 


40 


Differential equations form a vast field of study, and only the very 
beginnings have been described in this chapter. Their importance lies in 
the frequency with which they crop up as central relationships within 
mathematical models. You will see much more of differential equations if 
you undertake further studies in mathematics. 


Equations of the form dy/dx = f(x) and dy/dx = f(x)g(y) do not exhaust 
the types of first-order differential equations for which analytical solutions 
can be found. For example, the differential equation 


dy 
= 4.1 


whose direction field you studied with the computer in Section 4, has the 
general solution 


y = Ae* —x—1, 


where A is an arbitrary constant. Higher-order equations can also be 
solved analytically. For example, the second-order differential equation 
d*y ; 
—— + yy y = 0, 
ques 


where w is a non-zero constant, has the general solution 
y = Acos(wt) + Bsin(ut), 


where A and B are arbitrary constants. Note the presence of two arbitrary 
constants here, for the general solution of a second-order equation. To 
obtain a particular solution in this case, two initial conditions are required, 
perhaps corresponding to a known position and velocity at a given time. 
This is typical behaviour; the number of arbitrary constants that feature in 
the general solution is the same as the order of the differential equation. 


Where analytical solutions cannot be found, numerical methods may again 
be employed, for any order of differential equation. Euler’s method is just 
the first in a long line of methods which can be used for this purpose. 


We have so far considered equations that involve derivatives of a function 
of one independent variable, say y = f(x). It is also possible to develop 
calculus for functions of two or more independent variables, say z = f(z, y). 
Here there are two types of first-order derivative of z which may be taken: 
with respect to x, denoted 0z/0zx, or with respect to y, denoted 0z/Ody. 
These are called partial derivatives, in contrast to the ordinary derivatives 
considered in this block. Equations relating such derivatives, which again 
crop up within significant mathematical models, are called partial 
differential equations. As with ordinary differential equations, these may 
be of first or higher order. While some can be solved analytically, it is 
often the case that numerical methods are required to find their solutions. 


As a final remark, it is possible to combine these ideas with those of 
vectors, which you saw in Chapter B3, and hence to arrive at the topic of 
vector calculus. Both ordinary and partial differential equations feature 
here. 


Summary of Chapter C3 


In this chapter you were introduced to first-order differential equations. 
You saw how to solve equations of the form dy/dx = f(x) by direct 
integration, and equations of the form dy/dz = f(x)g(y) by separation of 
variables. The more general type of equation dy/dx = f(x,y) can be 
solved numerically. 


The continuous exponential model, embodied in the differential equation 
dy/dx = Ky (where K is a constant), has numerous applications. You saw 
it applied to radioactive decay and population change. If data are available 
in such a case, then values may be found for the parameters (such as K) 
from a log-linear plot. 


Learning outcomes 


You have been working towards the following learning outcomes. 


Terms to know and use 


Differential equation, order, first-order, solution, general solution, 
particular solution, initial condition, initial values, initial-value problem, 
direct integration, implicit solution, explicit solution, implicit 
differentiation, separable, separation of variables, exponential 
growth/decay, decay constant, half-life, proportionate growth rate, 
doubling time, discrete/continuous exponential model, log-linear plot, 
direction field, Euler’s method, step size. 


Notation to know and use 


y(a) = b (to express an initial condition on the variable y). 


Mathematical skills 


© Be able to recognise differential equations which can be solved by 
direct integration, and to solve suitable equations by this method. 


© Be able to recognise differential equations which can be solved by 
separation of variables, and to solve suitable equations by this method. 


© Knowing the general solution of a differential equation, be able to use 
an initial condition to find a particular solution. 


© Know the importance of checking a solution of a differential equation 
by substituting it back into the equation. 


© Be able to write down the general solution of a differential equation of 
the form dy/dx = Ky, where K is a constant. 


© Given data for the variables in such a differential equation, use a 
log-linear plot to estimate values for the parameters (such as K). 


© In an exponential decay process, be able to find the half-life from the 
decay constant, or vice versa. 


Modelling skills 


© Interpret predictions from an exponential model in real-world terms. 


Al 


Summary of Block C 


This block has introduced calculus, which is the cornerstone of 
mathematical studies at this level. Facility with calculus comes at the 
price of much exposure to it and practice with it, so do not be perturbed if 
you feel that you have yet to grasp fully some of the important points. 


Chapter Cl dealt with the topic of differentiation, starting with an 
informal explanation in terms of rates of change and the gradients of 
graphs, then providing a formal definition involving a limit. The 
consequences of this definition included a table of standard derivatives and 
various rules for differentiating sums and constant multiples of functions, 
products of functions, quotients of functions and composite functions. The 
alternative Leibniz notation was introduced. Differentiation was applied to 
the process of sketching graphs and to optimisation. 


Chapter C2 started from the basis that integration is both the reverse 
process to differentiation and a process in its own right. A table of 
standard indefinite integrals was drawn up, corresponding to that for 
derivatives, and rules were deduced for integrating sums and constant 
multiples of functions. The role of basic algebraic manipulation was 
stressed, as a means of modifying the expression for a function until it is in 
a form where its integral can be recognised, from the table or otherwise. 
Integration was applied to analyse one-dimensional motion, and in 
particular motion with constant acceleration, for which specific formulas 
were derived connecting time, position and velocity. It was shown that 
definite integration provides a method for finding areas under graphs. 


Chapter C3 involved in significant part revision and reinforcement of the 
material in the previous two chapters. Integration is central to the two 
strategies put forward for solving differential equations, and differentiation 
is needed in order to check solutions. (Also, in the shape of implicit 
differentiation with the Chain Rule, it permitted an explanation of why 
the separation of variables method works.) The principal modelling 
application was to exponential growth and decay, though other contexts 
were introduced also. 


This block has been but a brief initial tour of the possibilities provided by 
calculus, and you can expect to come across it again on a frequent basis if 
you continue with mathematical studies. This will provide reinforcement 
to your understanding of the basic ideas and results, as well as further 
appreciation of its use in modelling. 


Solutions to Activities — : 


Solution 1.1 


(a) 


The differential equation is dy/dr = 5/(2,/2). 
The derivative of the given function, 
y = 5/e = 5xr1/?, is 


d te 


— (5gl/2) = = 
7, (ee )}= $2 Wr 


which is the same expression as the right-hand 
side of the differential equation, as required. 


The differential equation is 

dy/dx = 3e%* cos(e®*). In order to differentiate 
the given function, y = sin(e®”), we apply the 
Composite (Chain) Rule. The function can be 
written as 


y=sinu, where u=e>’, 
for which 

d d 

= =cosu and = = 3c", 
Hence we have 

dy _ dy du 

dx dudz 

= (cos u)(3e°*) 


= 3e** cos(e**), 


which is the same expression as the right-hand 
side of the differential equation, as required. 


The differential equation is 

dy/dx = y(1 — 2tan(2z)). In order to 
differentiate the given function, y = e® cos(2z), 
we apply the Product Rule. The function can be 
written as y = uv, where 


esc — and — » = cvest2z), 


for which 
du = dv 
a —e€E and es =— —2 sin(22). 


Hence we have 
dy du dv 
de dx dex 
= (e”)(cos(2x)) + (e”)(—2 sin(2z)) 
= e”(cos(2z) — 2sin(2z)). 
Now we substitute y = e” cos(2z) into the 
right-hand side of the differential equation, to 


obtain 
y(1 — 2tan(2zx)) = e” cos(2z)(1 — 2 tan(2z)) 


e*(cos(2x) — 2sin(2z)). 


This is the same expression as was obtained 
above for the derivative, as required. 


(d) The differential equation is dy/dx = —3(y +1)’. 
In order to differentiate the given function, 
y = (1-—<2)/(1+ 2), we apply the Quotient Rule. 
The function can be written as y = u/v, where 


“=n i= 7 “and 7=1+ 2, 


for which 
du dv 
da ae 
Hence we have 


dy 1 PP du 
deg de ade 
(1 + x)(—1) — (1 — 2)(1) 

(1 + x)? 

ose 

(1+)? 

Now we substitute y = (1 — x)/(1+ 2) into the 
right-hand side of the differential equation, to 
obtain 


This is the same expression as was obtained 
above for the derivative, as required. 


Solution 1.2 

Substitute the values x = 1 and y = 7 into the 

expression y = 5,/x+c. This gives 
7=5V1+c=5+¢e, 


so c = 2. The required particular solution is therefore 


y = 5/xr +2. 


Solution 1.3 
The direct integration method can be applied to a 


differential equation of the form 
d(dependent variable) 

d(independent variable) 

Equations (a), (c) and (e) are of this form, so direct 


integration can be applied to these. Equations (b), 
(d) and (f) are not of the appropriate form. 


Be 


= f (independent variable). 
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Solution 1.4 


(a) The general solution of dy/dx = e?” is 


y= fe dx 


a Se 
= Z€ + C, 
where c is an arbitrary constant. 


(Check: Differentiating y = $e?” + c gives 


dy/dx = e**, as required.) 


i 
2 


(b) The initial condition is y= 2 when z = 0. 
Putting x = 0 and y = 2 into the general 
solution, we obtain 


aie ee 
a= 5¢ +c= 4 +€. 


Hence c = S, and the required particular 
solution is 


—— = 
ee: — 


Solution 1.5 


The general solution of dx/dt = sint is 


x= [ sintat 


= — cost + ¢, 
where c is an arbitrary constant. 


(Check: Differentiating 7 = — cost + c gives 
dx/dt = sint, as required.) 


The initial condition is x = 1 when t = 0. Putting 
t = 0 and zx = 1 into the general solution, we obtain 


1=-—cosO0+c=-—1+e. 


Hence c = 2, and the solution of the initial-value 
problem is 


© = 2 — cost, 


Solution 1.6 
(a) The general solution of dm/dt = —8t is 


m= | (—8t) dt 


= —4p +<, 
where c is an arbitrary constant. (This general 


solution can be checked by differentiation.) 


(b) The initial condition is m = 200 when t = 0. 
Putting ¢t = 0 and m = 200 into the general 
solution, we obtain 


200 = 0+. 


Hence c = 200, and the required particular 
solution is 


m = 200 — 4t?. 
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(c) When m = 100, the corresponding value of t is 
given by 


100 = 200 — 4t?;_ that is, t? = 25. 


Hence t = 5 when m = 100. The rocket’s fuel is 
exhausted after 5 seconds. 


Solution 2.1 


(a) The given equation is of the form H(y) = F(z), 
with H(y) =y+siny and F(x) = x+e* —1. It 
follows from result (2.3) that 


d ite d 
~(y+siny)— = piled eF = 1), 


dy a ar 
which leads to 
dy 
1 —=1+e’. 
(1 + cosy) a +e 


Dividing through by 1 + cosy (which is non-zero 
for —7 < y < 7) gives 

7 322 

dx 1+ cosy’ 
as required. 


(b) At (0,0), the slope of the curve is 


dy 1+e° 2 
ae Se eee Se Sp 
dx 1+ cos0 2 


Solution 2.2 


The given equation is of the form H(y) = F(x), with 
H(y) =tany and F(z) =z. It follows from 
result (2.3) that 


which leads to 


d 
sec? y 2 =. 


dx 

Dividing through by sec? y (which is non-zero for 
—50 <=. $7), and using the trigonometric identity 
sec? y = 1+ tan” y, we have 

dy _ 1 
dy 1+tan*g 
However, tan y = x, and so 

ee 

dx 1+?’ 


as required. 


Solution 2.3 


(a) The differential equation is (1/y?)dy/dx = —1. 
Comparing with result (2.4), we have 
h(y) = 1/y? = y~? and f(x) = —1. Hence the 
general solution of the differential equation is 
given by | 


[view [vac 


that is, 
tee Ee, 


where c is an arbitrary constant. 


(b) The explicit solution is obtained by making y the 
subject of the above equation. In this case, that 
is achieved by taking the reciprocal of each side 
and multiplying through by —1. The outcome is 


1 


’ 
i eee © 


y= 


where c is an arbitrary constant. (Since y > 0, 
we must have x —c > 0. There is some 
discussion of such restrictions on variables later 
in the section.) 


This answer could equally well be written as 
1 
r+’ 


where c is an arbitrary constant (in which case, 
we must have x +c > 0). 


y= 


Solution 2.4 


Equations (a), (b), (c) and (e) are of the appropriate 
form, but equations (d) and (f) are not. 


(a) f(t)=2, gly) =cosy 
(b) f(t) = beast -olep= t 
(c) f(t)=1, g(p) =pcosp 


(ce) fy)=1lt+y, g(v)=v 

These are the most likely choices for f and g, but 
they are not the only correct possibilities. For 
equation (a), for example, you could have chosen 
f(x) = 2a and g(y) = 4 cosy. 


Equation (b) can also be solved using direct 
integration; this is the most direct method of 
solution. 


_ SOLUTIONS TO ACTIVITIES 


Solution 2.5 


(a) 


The differential equation dy/dx = —x/y has a 
right-hand side of the form f(x)g(y), where 
f(x) = —x and g(y) = 1/y. Dividing both sides 
of the differential equation by g(y) means 
multiplying both sides by y. This gives 
dy _ 

Je 

(Other choices for the functions f and g are 


possible, but all choices lead to the same general 
solution. ) 


From this we obtain, using result (2.4), 


J vay = | (x) ae. 


(The solution thus far can be shortened by using 
the shortcut described in the text before 
Activity 2.4 on page 20. In that case, the 
solution to this point would simply read: 


Separating the variables in dy/dx = —x/y, we 
obtain 


[vay = [(-2)ae,) 


Carrying out the two integrations, we have 
sy — in? + C, 


where c is an arbitrary constant. This implicit 
form of the general solution can be manipulated 
to make y the subject, giving 


we 4 2c — 27, 


or equivalently 


a= fh 25 


where b = 2c is another arbitrary constant. We 
take the positive square root for y in order to 
satisfy the restriction y > 0 given in the 
question. 


(This expression for the general solution could 
be checked by substituting it into the differential 
equation. ) 


The initial condition is y = 1 when z = 0. 
Putting z = 0 and y = 1 into the general 
solution, we have 


1= /b—02 = Vb. 


Hence b = 1, and the required particular solution 


is 
y= V1—-2?. 


(For this expression to be well-defined, the range 
of values of x must be restricted. We shall 
return to this problem later in the section.) 
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Solution 2.6 


(a) 


The shortcut is adopted here. The differential 
equation is dx/dt = tsecx, so we can choose 
f(t) =t and g(x) = seca. Separating the 
variables in the differential equation, we obtain 


1 
/ da = | tat 
sec x 


that is, 


[coszae = | tat. 


After integration, we have 


an = St +c, 


where c is an arbitrary constant. Since 
—in <—2< ST, making the dependent variable x 
the subject of the equation, we obtain 


x = arcsin($t* + c) 
for the general solution in explicit form. 


The initial condition is x = 0 when t = 0. 
Putting t = 0 and x = 0 into the general 
solution gives 


0 = arcsin(0 + c) = arcsinc. 
Hence c = 0, and the required particular solution 
is 


x = arcsin($t*). 


Solution 2.7 


(a) 


Separation of variables in the differential 
equation gives 


fostom ecb 


Hence we have 
ty? =-+6C, 
r 


or equivalently (since v > 0) 


k 
v= a(S +e), 
r 


where c is an arbitrary constant. This is the 
general solution of the differential equation. 
The initial condition is v(R) = vg. Putting 
r= Rand v = vp into the general solution, we 
have 


k 
VR = 2(5 +e): 


Hence c = $v; — k/R, and the required 
particular solution is 


1 7 
v= s/o +k (= - 5), 


(c) As r becomes large, 1/r — 0, so the velocity 


tends towards the limiting value 


| 2k 
v2, — : 
If vp > \/2k/R then v2 > 2k/R, and so 


Zk 
2 
t= = > 0. 
ae 
Therefore, if up > ,/2k/R, the probe continues 
to travel away from the Earth indefinitely, since 
its velocity is never zero. 


On the other hand, if ve < \/2k/R then 

v?, < 2k/R, and so no real value for the velocity 
is obtained at sufficiently large distances from 
the Earth. In this case the probe does not travel 
that far away; it attains a maximum distance 
from the Earth and then falls back. 


(An alternative approach to this question would 
be to examine where the velocity of the probe 
becomes zero. This occurs where 


which gives a positive value for r only if 


UR < 1/2k/R. Hence the velocity of the probe is 
never zero if vg > \/2k/R.) 


The critical launch velocity, vp = \/2k/R, is 
called the Earth’s escape velocity. Its 
significance is that, according to the model, a 
probe launched with this or a greater velocity 
will escape completely from the Earth, whereas 
a probe launched with a lesser velocity will 
eventually fall back towards the Earth. (The 
value of the Earth’s escape velocity is 
approximately 1.12 x 104¢ms~!, that is, about 
40000 kph or 25000 mph.) 


Solution 3.1 


(a) The differential equation is 


dm/dt = —km (m > 0). Separating the 
variables, we obtain 


| cam=- [rat 
m 


Since m > 0, this gives 
Inm = —kt+c, 


where c is an arbitrary constant. This is the 
general solution in implicit form. 


Taking the exponential of each side, we obtain 


——— e ktte 


ae efe Kt = Ac. 
where A = e€¢ is a positive but otherwise 


arbitrary constant. This is the general solution 
in explicit form. 


(b) The initial condition is m = mo when t = 0. 
Putting t = 0 and m = mp into the general 
solution, we have 


mo = Ae® =A, 
Hence A = mo, and the required particular 


solution is 


Tt = nce. 


Solution 3.2 


Using equation (3.5), the decay constant k 
corresponding to a half-life J’ = 2.62 hours is 


inZ- 2 
SS SS SS eS a h = ia ci 
F 5 6) 0.265 hour ~ (to 3 s.f.) 


Hence, from equation (3.4), the proportion of the 
original amount remaining after six hours is 
m(t + 6) 
m(t) 


k 


ee ee” 0.204 (fo 3 5-f). 


So about 20% of the original amount of silicon—31 
remains after six hours. 


Solution 3.3 


Using equation (3.5), the decay constant k 

corresponding to a half-life 7’ = 5570 years is 
ee In 2 2 In 2 
f= SS 


At the present time, 85% of the original amount mo 
of carbon—14 remains, so we have 


k = 1.24 x 10-* year" (to 3 s.f.). 


p25 2 Seer 
™o 


where t (in years) is the time since the animal died. 
Taking the logarithm of each side, we obtain 

—kt = |n(0.85); 
that is, 


In(0.85) 


k 


In(0.85) 
= eee a to 3s.f.). 
1.24 x 10-4 eee Set 


So the animal died approximately 1300 years ago. 


“SOLUTIONS TO ACTIVITIES — 


Solution 3.4 | 
(a) The completed table of data is below. 


0.73 0.56 0.42 0.29 0.24 
—0.31 —0.58 -0.87 —1.24 —1.438 


m/mo 


(b) The graph of In(m/mo) against ¢ for the data in 
part (a) is shown in Figure S.1. This figure also 
shows a straight line (chosen by eye) through 
the origin which fits the data points well. (This 
good fit confirms that the radioactive decay of 
uranium—239 can be modelled satisfactorily as 
an exponential decay process.) 


In(m/mo) 


—0.5 4 


—1.5H 


Figure S.1 


(c) The straight line passes through the origin and 
the point (50, —1.48), so its slope is 


1.48 


Hence the decay constant for uranium—239 is 
k = 0.0296 min~*. From equation (3.5), the 
corresponding half-life is 


T= In 2 = In 2 
k 0.0296 


= 23.4 min {te-3-s.£). 


AT 
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Solution 3.5 


(a) The initial population size is Py = 5 x 10°, and 
the proportionate growth rate is 
kK = 0.03 year~+. Hence, from equation (3:7), 
the population size P after time t years is 


P = Pye**t =5 x 1089-9, 
After 10 years the population size is 
P.=5.x:10"'¢"? = 6.75 x 10° {to 3 sf), 
whereas after 100 years, the population size is 
P=bhxl0’e = 100x440 Goa 3af). 


So, to three significant figures, the population 
size after 10 years is 6.75 million, and after 100 
years it is 100 million. 


(b) The population size reaches 50 x 10° at the time 
t years given by 


50 x 10° = 5 x 10%e 9%. that is, e%°%* = 10. 


Taking the logarithm of each side, we obtain 
0.03t = In 10; that is, 


t = ~°ln10 = 76.8 (to 3 sf.) 


Hence it takes almost 77 years for the 
population size to reach 50 million. 


Solution 3.6 


(a) From equation (3.8), the population will double 
inatime JT given by © 


9 — PUt+T) _ xr 
P(t) 
Taking the logarithm of each side, we obtain 
KT = 1n2, so the doubling time is 


In 2 
f ee ees 
K 


(b) For the population described in Activity 3.5, the 
proportionate growth rate is K = 0.03 year~!. 


Hence the doubling time is 


In 2 
‘ = = 23.1 years (to 3 s.f.). 


Solution 3.7 


(a) The table below gives the values of the time t 
(in years since 1950), the population size P (in 
millions) and In P, to three significant figures. 


398 442 555 689 851 1014 
InP | 5.88 6.09 632° 684 6.75 6.92 


The corresponding log-linear plot is shown in 
Figure 8.2. 
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Figure 8.2 


The data points are fitted reasonably well by the 
straight line drawn. This confirms that the 
exponential model is appropriate to describe the 
Indian population for the years 1950-2000. 


The intercept on the (In P)-axis is 5.89; this is 
the value of In Py for the model. Since 
In Po = 5.89, we have 


Fy =e" = Set fie 3 of). 


As well as the point (0,5.89), another point on 
the straight line is (50,6.94). The slope of the 
straight line is therefore 

6.94 — 5.89 

50 — 0 

So our estimate for the proportionate growth 
rate K is 0.021 year+, and for the initial 
population Po the estimate is 361 (million). 


= 0.021. 


Using the estimates obtained for K and Po, the 
exponential model for the Indian population 
(1950-2000) is 


P = 361e°-°?1%, 


The year 2020 corresponds to t = 70, so our 
estimate for the population in that year is 


P(70) = 36ie"°""*" = 1570 (46 9 92.). 


So our estimate for the population of India in 
2020 is 1570 million, that is, 1.57 billion. 


Solutions to Exercises 


Solution 1.1 
The differential equation is dy/dz = y?. The 


derivative of the given function, y = —1/z, is 
d = eee = eae | 
= (-2)--z¢ pt as 


Now we substitute y = —1/z into the right-hand side 
of the differential equation, to obtain 


2 
Pee 
xL x2 


This is the same expression as was obtained above 
for the derivative. Hence y = —1/z (x > 0) isa 
solution of the given differential equation. 


Solution 1.2 
(a) The general solution of dy/dx = ./x = x1/? is 


y= fo? dx 


me Se 
= gv / + C, 
where c is an arbitrary constant. 


(b) The initial condition is y = 5 when z = 4. 
Putting x = 4 and y = 5 into the general 
solution, we obtain 


ee 3/2 ae 
b= 2x 4? 4+c= Be, 


Hence c = —i, and the required particular 
solution is 
Sage fos gee | 


(c) When x = 25, we have 
y= 2x 259/22 
=2x125-4 

= 5(250 — 1) = 83. 


Solution 1.3 
(a) The general solution of du/dz = cos(2z) is 


wih / cis oy de 


= ;sin(2r) +c, 


where c is an arbitrary constant. Using the 


initial condition, u = —2 when x = aT, we 
obtain 

—2= $sin($7) +c=$+4¢. 
Hence c = —2, and the solution of the 


initial-value problem is 


u = 5sin(2x) — 2. 


(b) The general solution of dx/dt = 1/t (t > 0) is 


1 
c= [iat 
t 
=Int+c, 


where c is an arbitrary constant. Using the 
initial condition, x = 3 when t = 1, we obtain 


o> ini+-c=c. 
Hence c = 3, and the solution of the initial-value 
problem is 

x=Int+3. 


Solution 1.4 | 


(a) The general solution of the differential equation 
dr/dt = —k is 


r=— kat 


= —kt+c, 
where c is an arbitrary constant. 


(b) Since the initial radius of the snowball is 100 cm, 
the initial condition is r = 100 when ¢ = 0. 
Substituting this into the general solution, we 
obtain c = 100. Hence the required particular 
solution is 


r = 100 — kt. 


(c) The initial volume of the snowball (in cm?*) is 
a x 100°. If we denote the radius of the 
snowball after 2 days by R, then 


-aR° — 5 x <7100°, 
from which we find 
R=27'/3 x 100. 


Using the condition r = R = 2~1/% x 100 at 
= 2 in the particular solution of the differential 
equation, we obtain 


2-1/3 x 100 = 100 — 2k. 
It follows that 
k = 50(1 — 271/3) ~ 10.315. 


(d) The snowball disappears when r = 0; that is, 
when 100 — kt = 0. This occurs when 


= 100 _ 2 

SS ee 
It therefore takes about 9.7 days for the 
snowball to disappear. 


= 9.695. 
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Solution 2.1 


(a) 


(b) 


The given equation is of the form H(y) = F(z), 
with H(y) = y+ lny, where y > 0, and 

F(x) = (1 — 2z)e°*. It follows from result (2.3) 
that 

d dy d 3 

— Iny)— = —((1 — 2x)e”* 

gy + Iny) Se = © ((1— 22)e%), 

which (using the Product Rule on the right-hand 
side) leads to 


(: + ~) ay = (—2)(e°") + (1 — 22)(3e**); 


y/ dx 
that is, 
ld 
YS 4 ee 
ya 


Multiplying through by y/(y + 1) gives 


dy  y(1—6z)e** 
a 


) 


as required. 
At (0,1), the slope of the curve is 
dy — {i= 2 


de --- 442 


Solution 2.2 


(a) 


The differential equation is u? du/dx = x*. 


Comparing with result (2.4), we have h(u) = u? 
and f(x) = x*. Hence the general solution of the 
differential equation is given by 


[tana f xtae; 


where c is an arbitrary constant. Solving this 
equation for wu gives the explicit form of the 
general solution as 


1/3 
us (20° +0)", 
where b = 3c is an arbitrary constant. 


The initial condition is u = 1 when zx = 1. 
Putting x = 1 and u = 1 into the general 
solution, we obtain 


1= (+9) 


Hence 6 = 2, and the required particular 
solution is 


1/3 


1/3 
w= (B24 2)"" 


Solution 2.3 


(a) 


The differential equation is 


dy/dx = —y° (y > 0). Separating the variables, 
we obtain 


[(-p)a- fue 


oe 
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that is, 


—~ =F +6, 
2y? 
where c is an arbitrary constant. To make y the 
subject of this equation, first take the 
reciprocals of both sides and then divide through 
by 2, to obtain 
——— 
2(r@+c) 
Since y > 0, the general solution in explicit form 
is therefore 
1 


ere | 


The initial condition is y = 5 when x = 0. 
Putting + =C and 7 = 5 into the general 
solution, we obtain 


1 1 


—=— = 
Hence c = 2, and the required particular solution 
is 
= 1 
~ aera) 


The differential equation is dx/dt = tx (x > 0). 
Separating the variables, we obtain 


| cae [eat 
g 


that is, 
= St? + ¢, 


where c is an arbitrary constant. To make zx the 
subject of this equation, recall that 

exp(Inz) = x. Hence, taking exponentials of 
both sides, we have 


sa ee 


2 2 
= ee" 4 = Ae’ /?; 


where A = e°. The constant A must be positive 
(since e° > 0 for any choice of c) but is otherwise 
arbitrary. The general solution for x > 0, in 
explicit form, is therefore 


2 
z= Ae’ sa 


where A is a positive but otherwise arbitrary 
constant. 


The initial condition is x = 2 when t = 0. 
Putting t = 0 and x = 2 into the general 
solution, we have 


2 =Ae”. 
Hence A = 2, and the required particular 


solution is 


= det /2. 


Solution 2.4 


The volume V is decreasing, and the rate of change 
of volume dV /dt is proportional to the surface 
area A. In other words, 


dV 
—=-—kA 
dt 


where k is the (positive) constant of proportionality. 
Putting V and A in terms of r, this differential 
equation becomes 


s (4rr*) = —Arkr?. 


Hence, using the Chain Rule, we have 
d 
Anr? — = —Arkr?, 
dt 


which gives 
dr _ 
a 
Solution 2.5 


(a) The differential equation is 
dh/dt = —kh‘/? (h > 0). Separating the 
variables, we obtain 


pe dh = [-» dt: 


that is, 
2ni/? = —kt +c, 


where c is an arbitrary constant. Hence the 
explicit form of the general solution is 


h = 4(c— kt)’. 


(b) The initial condition is h = ho when t = 0. 
Putting t = 0 and h = hg into the general 
solution, we obtain 


ho = =c. 


Now the equation 2h!/? = —kt + c above shows 
that c must be positive, since we have il > 0, 
k >0Oandt>0. Hence c= ah. and the 
required particular solution is 


k= Ce = kt) 
(c) The tank empties to the level of the hole when 
h = 0. This occurs when 
he! a Skt =: 
that is, at time 


_ oni? 


t 
k; 


SOLUTIONS TO EXERCISES © 


Solution 3.1 


(a) 


The credit balance of the account increases at a 
proportionate rate 0.05, so the balance c at 
time t satisfies the differential equation 

dc 
By comparison with result (3.9), the general 
solution is 

= Aa 


where A is an arbitrary constant. (However, we 
require A > 0 for the purposes of the model.) 


If the initial balance is £100, then the initial 
condition is c = 100 when t = 0. Substituting 
these values into the general solution, we obtain 
100 = Ae° = A, 
so the required particular solution is 
c= 100e°*. 
The credit balance after one year is therefore 
el = 1002 = 105.13, 


to the nearest penny. At the end of one year, the 
balance of the account is £105.13. 


In this part we use the formula given in the 
margin note, with cp = 100. 


The account balances (in £) after one year, for 
an interest rate of 5% per year (r = 0.05) 
compounded annually (n = 1), six-monthly 

(n = 2) and quarterly (n = 4), are respectively 
as follows: 


100 x 1.05 = 105; 

100 x (1.025)” = 105.06; 

100 x (1.0125)* = 105.09. 
The interest accrued on £100 after one year at 
5% interest is therefore £5 (compounded 
annually), £5.06 (compounded six-monthly), 


£5.09 (compounded quarterly) or £5.13 
(compounded continuously). 
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analytical solution 38 


carbon dating 28 
continuous exponential model 32 


decay constant 26 

estimation 29 
differential equation 

first-order 6 

separable 19 
direct integration 11, 12 
direction field 38 
discrete exponential model 32 
doubling time 32 


Euler’s method 39 
explicit solution 15 
exponential model 32 


function of two variables 38 
general solution of a differential equation 9 
half-life 28 


implicit differentiation 16 
implicit solution 15 
initial condition 10 
initial value 10 
initial-value problem 10 


log-linear plot 33 
order of a differential equation 6 


particular solution of a differential equation 9 
proportionate growth rate 31 


separation of variables 19 

solution curve 13 

solution of a differential equation 6 
step size 39 
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